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ABSTRACT 

Observations of the Sun with the LOWL instrument provide a homogeneous set of solar 
p-mode frequencies from low to intermediate degree which allows one to determine the 
structure of much of the solar interior avoiding systematic errors that are introduced 
when different data sets are combined, i.e., principally the effects of solar cycle changes 
on the frequencies. Unfortunately, the LOWL data set contains very few of the lowest- 
degree modes, which are essential for determining reliably the structure of the solar 
core - in addition, these lowest-degree data have fairly large associated uncertainties. 
However, observations made by the Birmingham Solar-Oscillations Network (BiSON) 
in integrated sunlight provide high-accuracy measurements of a large number of low- 
degree modes. In this paper we demonstrate that the low-degree mode set of the LOWL 
data can be successfully combined with the more accurate BiSON data, provided the 
observations are contemporaneous for those frequencies where the solar-cycle-induced 
effects are important. We show that this leads to a factor-of-two decrease in the error 
on the inferred sound speed in the solar core. We find that the solar sound speed is 
higher than in solar models for r < O.2R0. The density of the solar core is, however, 
lower than that in solar models. 
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1 INTRODUCTION 

Observations of the Sun, covering many years, have now 
provided us with accurate measurements of solar p-mode 
frequencies which impose severe constraints on the structure 
of the Sun. 

The LOWL instrument (LOWL is an abbreviation for 
low degree with degree denoted by L) has, for the first time, 
provided us with a uniform set of frequencies from low to 
intermediate degree I. This has allowed detailed inversions 
for the structure of much of the Sun's interior (cf., Basu 
et al. 1995; 1996a; 1996b) to be performed. The available 
data from the instrument now span a period of more than 1 
year, and most of the mode- frequency determinations have 
relative errors as small as a few parts in 10®. 

The LOWL data set contains very few of the lowest- 
degree modes, i.e., for < / < 2, and those that are present 
have much larger errors than their higher-/ counterparts. 
Since only the low-degree modes penetrate to the solar core, 



large uncertainties in the inferred structure of the solar core 
arise from inversions which rely solely on these data. Ob- 
servations made by the Birmingham Solar-Oscillations Net- 
work (BiSON) do, however, provide very accurate measure- 
ments of low-degree modes (cf. Elsworth et al. 1994). These 
modes have to be combined with intermediate and high- 
degree modes from other sources before they can be used to 
infer the solar structure by inverting the observed frequen- 
cies - it would therefore seem logical to combine the BiSON 
and LOWL solar data in order to provide a more reliable 
determination of the structure of the solar core. 

Great care must be taken in combining data from dif- 
ferent sources. Systematic differences may arise from differ- 
ences in instrument characteristics and analysis techniques, 
or from temporal variations in the solar p-modes if the differ- 
ent observations are not contemporaneous. Such differences 
may be interpreted by an inversion as being solar in ori- 
gin, and in particular in the present context as arising from 
spatial variations in the solar interior. 
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Table 1. Summary of data sets 
Name Summary details 

1 BiSON-8 8-month BiSON spectrum; maximum 

likelihood analysis 

2 BiSON-G Five 2-month BiSON spectra, Gaussian- 

analysed 

3 Lower frequency data from 32-month BiSON 
spectrum; maximum likelihood analysis 

4 Comparison BiSON data taken at high solar 
activity; Gaussian-analysed 

5 LOWL 1-ycar spatially-resolved LOWL data 

6 Best Set Combination of BiSON sets 1-3 and LOWL 

The issue of combining non-contemporaneous data is 
particularly important, because it is known that the p- 
mode frequencies change along the solar cycle (Libbrecht & 
Woodard 1990; Elsworth et al. 1994). The difficuhies in de- 
termining the solar core structure from inhomogeneous data 
sets have been discussed by Gough & Kosovichev (1993), 
Gough, Kosovichev & Toutain (1995) and Basu et al. (1995, 
1996a). The solar-cycle variations are believed to be the re- 
sult primarily of variations in the surface properties of the 
Sun. It may therefore be possible that the temporal vari- 
ations can be removed in the inversion in the same way 
as other surface uncertainties (cf. Dziembowski et al. 1991; 
Kosovichev et al. 1992) but, as was pointed out by Basu et 
al. (1996a), there remain significant problems in the combi- 
nation procedure and the resulting inversions may be mis- 
leading. 

We have largely avoided such problems in the present 
work by combining near-contemporaneous data from the Bi- 
SON network and the LOWL instrument, taken at a low- 
activity phase of the solar cycle. We also use some low- 
frequency BiSON data based on observations over several 
years, but such low- frequency modes will be relatively in- 
sensitive to solar-cycle variations in the near-surface layers 
(cf. Libbrecht & Woodard 1990). We also present results 
of combining LOWL with BiSON data collected at times 
of high solar activity, to illustrate the dangers inherent in 
inverting non-contemporaneous data. 



2 THE DATA 

The data used for the analysis come from two sources: the 
LOWL instrument based on Mauna Loa (Tomczyk et al. 
1995) in the Hawaiian Islands; and the global, 6-station Bi- 
SON network (Chaplin et al. 1996a). Both systems employ 
the same basic physical principles - the use of an atomic 
standard - to measure the Doppler velocity shift of a so- 
lar Fraunhofer line formed by potassium atoms in the near 
infrared (770 nm). The LOWL instrument spatially resolves 
the visible solar disc, and is therefore sensitive to oscillation 
modes of degree up to / = 99. The BiSON instruments view 
the unresolved Sun, and are sensitive to modes of < / < 4. 
Although limited to the lowest-Z modes, the BiSON tech- 
nique is very stable and provides some of the highest-quality 
measures of those modes available. The combination of these 
two data sets should be very powerful given that: they result 
from observations made on the same line in the solar atmo- 
sphere; each are characterized by high intrinsic accuracies; 



and both were collected at the same epoch. 

The LOWL time series used for the analysis covers a 12 
month period beginning 1994 February 26 and has a tempo- 
ral duty cycle of 22 per cent. Further details of the LOWL 
experiment are given by Tomczyk et al. (1995). The data 
were analysed up to 3.5 mHz, which thus set the upper limit 
to the frequency range used in this paper, while the lower 
limit was set by the visibility of the modes in the data set. 
The fitting procedure, which is a maximum-likelihood tech- 
nique applied to the complex Fourier transform of the time 
series, is described in Appendix 5 of Schou (1992). 

Several BiSON data sets have been used: these are sum- 
marised briefiy in Table 1. For reasons of convenience, we 
have used an 8-month BiSON spectrum within the time span 
of the LOWL set (henceforth referred to as BiSON-8), gener- 
ated from data collected between 1994 January 26 and 1994 
September 22. The duty cycle of the BiSON-8 time series 
was 80 per cent. The modes in the Fourier spectra gener- 
ated from the BiSON-8 time series were fitted by minimizing 
a maximum-likelihood function that assumed an underlying 
X distribution with two degrees of freedom in the frequency 
domain (e.g., Anderson, Duvall & Jefferies 1990; Toutain & 
Appourchaux 1994). 

In addition, we have fitted frequencies from five 2- 
month BiSON spectra (which are contemporaneous with the 
LOWL data). The duty cycle for the complete calendar year 
1994 was 78 per cent. The modes in the 2-month BiSON 
spectra were first smoothed before fitting to yield pseudo- 
Gaussian statistics in the frequency domain - these data are 
therefore referred to as the BiSON-G set. 

To get reliable determinations of the lowest-degree 
modes at frequencies below 1.8 mHz one requires spectra 
generated from more than 1 year's data, and for this rea- 
son we used results derived from the not-entirely contempo- 
raneous 32-month BiSON spectrum. These modes are not 
expected to be influenced significantly by the solar cycle, 
and we do in fact demonstrate that their inclusion in the 
analysis does not bias the results. Mode frequencies were 
obtained from a single, 32-month BiSON spectrum of obser- 
vations between 1992 October and 1995 June (Chaplin et 
al. 1996b). The fractional fill of useful data was 72 per cent. 
The spectrum was analysed in the same way as the 8-month 
spectrum. 

Finally, to show that variations in p-mode frequencies 
over the solar-cycle do have a deleterious effect on structural 
inversions, we have used a set of frequencies derived from 
BiSON observations made at times of high solar activity. 
These data are described by Elsworth et al. (1994). 

The fitted BiSON frequencies at Z = and 1 were, in 
general, more accurate than their LOWL counterparts; for 
1 = 2, the quality of the data were similar; while for 1 = 3, the 
LOWL frequencies were, by and large, better determined. 
We have derived a "Best Set" of frequencies by merging the 
available data, giving due and careful consideration to the 
comparative quality of the frequency determinations in each 
set (see later) . The LOWL and BiSON data sets are plotted 
in Fig. 1 in the form of an l-v diagram. Because the lower 
turning points of the modes in the combined data set span 
the range of most radii from the Sun's surface to its centre, 
and because the errors on the frequencies are very small, it 
is feasible to perform a meaningful inversion for the solar 
structure throughout most of the interior. 
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Figure 1. The l-u diagram of the data used. Panel (a) shows the entire range with lOOOo" error bars, while panel (b) shows just the 
low-degree modes. The triangles (dots in panel a) are LOWL modes; eircles are BiSON-8 modes; asterislts are 32-month BiSON modes. 



3 INVERSION TECHNIQUES AND SOLAR 
MODELS 

Solar oscillations can be described throughout most of the 
solar interior by equations describing linear adiabatic oscil- 
lations (cf. Unno et al. 1989). Chandrasekhar (1964) showed 
that these equations, along with appropriate boundary con- 
ditions, constitute a self-adjoint eigenvalue problem, which 
leads to a variational principle connecting the eigenfrequen- 
cies to the basic equilibrium state of the Sun. 

Our inversion for solar structure is based on linearizing 
the equations of stellar oscillation around a known reference 
model (cf., Dziembowski, Pamyatnykh & Sienkiewicz 1990; 
Dappen et al. 1991; Antia & Basu 1994; Dziembowski et 
al. 1994) . The differences between the structure of the Sun 
and the reference model are then related to the differences 
in the frequencies of the Sun and the model by kernels. 

Non-adiabatic effects give rise to frequency shifts (Cox 
& Kidman 1984; Balmforth 1992) which are not accounted 
for by the variational principle. Frequency shifts are also in- 
troduced by errors in modelling the underlying solar model, 
e.g., the effects of turbulent convection. Most of these un- 
certainties affect the surface layers of the models. In the 
absence of any reliable formulation, these effects have been 
taken into account in an ad hoc manner by including an ar- 
bitrary function of frequency in the variational formulation 
(Dziembowski et al. 1990). This can be justified because 
in the surface layers, the eigenfunction is largely indepen- 
dent of the degree I of the mode; thus, for spherically sym- 
metric perturbations, the frequency difference resulting from 
the surface effects, weighted by the inertia of the mode, is 
roughly a function of frequency only (Christensen-Dalsgaard 
& Berthomieu 1991). We thus express the fractional change 
in frequency of a mode in terms of fractional changes in the 
model parameters and also a surface effect. 

When the oscillation equation is linearised - under 
the assumption of hydrostatic equilibrium - the fractional 
change in the frequency can be related to the fractional 



changes in two of the model parameters. Thus, 



= / ^1.2 r dr + 
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+ 
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(1) 



(cf. Dziembowski et al. 1990). Here SuJi is the difference in 
the frequency uji of the ith mode between the solar data 
and a reference model. The functions /i and /2 are an ap- 
propriate pair of model parameters. The kernels KI2 and 
A'2,1 are known functions of the reference model which re- 
late the changes in frequency to the changes in /i and /2 
respectively; and Qi is essentially the inertia of the mode 
(Christensen-Dalsgaard 1986). The term in Fsut! results 
from the near-surface errors. 

The pair of variables (/i, /2) can involve several combi- 
nations of model parameters. As discussed below, our goal 
in the inversion is to isolate one of the variables while en- 
deavouring to ensure that our results are insensitive to the 
other variable. In this work we generally use {c^,p), c being 
adiabatic sound speed and p density, to invert for c^. To in- 
vert for density, one good pair to use is (p, Fi), where Fi is 
the first adiabatic exponent. Alternatively, one may assume 
the equation of state to be known and transform the depen- 
dence of the oscillations on Fi into a dependence instead on 
Y (Y being the abundance by mass of helium), to carry out 
inversion for the pair {p,Y). (As formulated, this also as- 
sumes that the heavy element abundances are known.) We 
have used this formulation, which has the advantage that 
the frequency dependence on Y is largely confined to the 
helium ionization zones. This is easier to suppress in the 
inversion for density than the dependence on Fi would be 
in an inversion using (p, Fi). Consequently we can achieve 
better resolution and/or error properties in our density in- 
version. However, this advantage is bought by assuming the 
equation of state to be known, a point to which we return 
in Section 5. 
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A number of different inversion techniques can be used 
for inverting the constraints given in equation (1). We have 
used the subtractive optimally locahsed averages (SOLA) 
method of Pijpers & Thompson (1992). 

The principle of the inversion technique is to form linear 
combinations of equations (1) with weights Ciiro) chosen 
such as to obtain an average of 5fi/ fi localized near r = 
ro while suppressing the contributions from 5/2 //2 and the 
near-surface errors. In addition, the statistical errors in the 
combination must be constrained. If successful, the result 
may be expressed as 

/K;(ro,r)^dr.5:c,(n,)^, (2) 

where the /C(ro, r), the averaging kernel at r = ro, is defined 



Table 2. Properties of solar models 
Model {Z/X)s Yc Yc T,_ 



pc r-d/R© Age 
10'' K g cm-3 Gyr 



JC{ro,r) = ^Ci{ro)Kl2{r) , 



(3) 



of unit integral, and determines the extent to which we have 
achieved a localized measure of Sfi/fi. In particular, the 
width in r of /C(ro,r), here calculated as the distance be- 
tween the first and third quartile point, provides a measure 
of the resolution^ . The precise implementation of the method 
is described by Basu et al. (1995, 1996c). 

A drawback of the SOLA method is the implicit as- 
sumption that the frequency differences and the errors as- 
sociated with them are correct. However, sometimes the ob- 
servational errors can be either under- or overestimated (e.g. 
because the modes have much lower line widths than modes 
adjacent in frequency), in which case the inversion result can 
be misleading. We have found that a reliable way to detect 
such modes is to do a regularized least squares (RLS) inver- 
sion first. We rejected modes which had more than a 3. So- 
residual after the fit. The SOLA inversion was performed 
on the weeded mode set. On examining some of the low- 
degree modes that had been rejected by the RLS procedure, 
we found that the fit to the mode concerned was usually 
poor. Thus the RLS is a reasonably reliable way of weeding 
out uncertain modes. There are some interesting instances 
where, even in contemporaneous data, the appearance of a 
particular mode was very different in each data set, presum- 
ably because of beating with noise. In addition, different 
analysis techniques can give rise to fitted frequencies which 
differ by more than the formal error. The details of the RLS 
implementation can be found in Antia & Basu (1994). 

Our inversion requires the use of a reference solar model. 
The model we use (Model S of Christensen-Dalsgaard et al. 
1996) was constructed with the Livermore (OPAL) equation 
of state (Rogers, Swenson & Iglesias, 1996). For tempera- 
tures higher than 10* K, OPAL opacities are used (Iglesias, 
Rogers & Wilson 1992), whereas at lower temperatures opac- 
ities from the tables of Kurucz (1991) are taken. The model 
incorporates the diffusion of helium and heavy elements be- 
low the convection zone. The surface heavy element ratio is 
Z/X = 0.0245 (Grevesse & Noels 1993). The model has an 



t The radii ri, r2, r^ of the first, second and third quartile points 
are defined such that J{ri) = 0.25, J{r2) = 0.5 and Jirs) = 
0.75, where J^{r) = I ICdr (recall that JC is unimodular). In 
practice this uniquely defines ri, r2 and rs for our kernels. 



Reference 0.0245 0.2447 0.6444 15.67 154.2 0.7115 4.6 
Test 0.0245 0.2457 0.6402 15.64 152.9 0.7124 4.52 

Notes: {Z/X)b is the present surface heavy element abundance 
ratio; Ye and Yc are the current envelope and central helium 
abundances respectively; Tc is the central temperature, pc the 
central density, and t^/Kq the radial location of the base of the 
convective envelope in the models used in this paper. 

age of 4.6 Gyr. We have also used another model as a proxy 
Sun to illustrate certain effects. The physical assumptions 
in the model are identical to those in the reference model, 
but it has a lower age of 4.52 Gyr. Some of the properties 
of the models are listed in Table 2. 

4 RESULTS 

4.1 Consistency of the frequencies 

To illustrate the effects of the solar cycle on the observed 
frequencies, the frequency differences between the LOWL 
and the BiSON-8 and BiSON high-activity data have been 
plotted in Fig. 2. The high-activity BiSON data show a 
systematic, frequency-dependent difference with respect to 
the LOWL data. The BiSON-8 data show no such trend. 
Since slowly varying, frequency- dependent terms in the dif- 
ferences are caused by perturbations near the solar surface 
(e.g. Gough 1990), this implies that while the BiSON-8 
and LOWL data show similar surface effects, the BiSON 
high-activity data do not. The differences between the high- 
activity BiSON and LOWL data are indicative of the in- 
fluence of the solar cycle. Consequently, if the low-degree 
LOWL data were to be replaced by the high-activity Bi- 
SON data, the combined mode set would have two different 
surface terms. 

This effect is seen quite clearly in Fig. 3(a), which shows 
differences between the LOWL, BiSON-8 and BiSON high- 
activity data and frequencies of the reference solar model. 
The frequency differences have been weighted with the cor- 
responding mode inertia in order to compensate for the 
fact that, for given a discrepancy in the model or mode 
physics, deeply penetrating modes are perturbed less than 
modes which have a shallow turning point (cf. Christensen- 
Dalsgaard & Berthomieu 1991). To show the general trend 
of the LOWL data, we have also plotted the frequency differ- 
ences of higher-degree modes. While the BiSON-8, LOWL, 
and 32-month BiSON data follow the same trend, the high- 
activity BiSON data lie above this. 

Owing to the manner in which the inversion techniques 
are implemented, they would be unable to isolate the two 
separate frequency dependent trends resulting from the com- 
bination of the high- and low-activity data. This would, of 
course, give rise to misleading inversion results (cf. Basu et 
al. 1995, 1996a). The BiSON-8 and LOWL data show no 
systematic differences, however, and we should therefore be 
able to combine these two data sets. 

Asymptotically, the frequency differences can be written 
as (e.g. Christensen-Dalsgaard, Gough & Perez Hernandez 
1988) 



S(w)5oj/u; = Hi{w) + H2(io), 



(4) 
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Figure 2. (a) The frequency diflferences between the high-activity 
BiSON and LOWL data sets, (b) The frequency differences be- 
tween the BiSON-8 and LOWL data sets. The error bars on each 
plot were generated by adding the respective BiSON and LOWL 
uncertainties in quadrature. 



where w — uj/{l + 0.5), and S{w) is a known function of the 
reference model. The function Hi (w) depends on the sound- 
speed difference between the Sun and the reference model 
and H2{ijj) is determined by differences at the surface. In 
Fig. 3(b) we have plotted the residuals after the functions 
H\ and H2 - each obtained by fitting the LOWL data - 
have been removed from the frequency differences. Note that 
the residuals for the high-activity BiSON data are much 
larger compared with those for the other sets, and follow a 
completely different trend. 

4.2 Consistency of inversion results 

Once the consistency of the different data sets had been 
established, we proceeded to invert for the sound-speed dif- 
ference between the Sun and the reference model, comparing 
results of different combinations of LOWL and BiSON data. 
The results of the inversions are shown in Fig. 4. To provide 
a common reference, all panels include the results of invert- 
ing the LOWL data only. 

First, in Fig. 4(a) we replaced the / = 0, 1, and 2 modes 
of the LOWL set with those of BiSON-8 or high-activity Bi- 
SON data. Only those modes which were common in all sets 
were used - this ensured that the resultant resolution of the 
inversions would be comparable. Since the mode sets used 
for the inversions were identical for I > 2, any differences in 
the results should be confined to the core region - this is in- 
deed found to be the case. As anticipated, the results of the 
inversion which utilized the high-activity BiSON set were 
markedly different from those based on data from the same 
epoch. This effect was first noted by Basu et al. (1996a), 
and is a result of the different surface effects in the high and 
low-activity data. As can be seen in the figure, the inversion 
results from the LOWL set and the LOWL plus BISON-8 
set are fairly consistent - the contemporaneous BiSON-8 and 
LOWL sets can therefore be satisfactorily combined. Simi- 
larly, inversion results [Fig. 4(b)] generated by substituting 
the low-degree data from the BiSON-G set show no basic 
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Figure 3. (a) The frequency difference between the different data 
sets and the reference model, in the sense (Sun) - (Model), scaled 
by the mode inertia Q„;, which is essentially the function S (cf. 
eq. 4) normalised by the acoustic radius of the star, (b) The resid- 
uals after the function Hi (ui) + H2 {ijj) has been removed from the 
scaled frequency differences. The tiny crosses are LOWL data for 
degrees 3 to 10; the filled triangles are LOWL data for degrees 
0—2; the empty circles are BiSON-8 data; the filled circles are the 
32-month BiSON data; and the filled squares are BiSON data 
obtained at high solar activity. 



inconsistency; thus, although fitted by a different technique, 
these mode frequencies can be used in a combined data set. 
As a final step, we supplemented the input data with 
those low-degree modes - from the BiSON-8 and 32-month 
BiSON sets - not present in the LOWL set. The results are 
shown in Fig. 4(c). Adding the extra modes does not change 
the general trend of the result, although there are noticeable 
improvements in resolution and errors. 
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results from a set where LOWL data for < Z < 2 have been 
replaced by BiSON-8 data; additional mode frequencies from the 
BiSON-8 and 32-month BiSON sets - not available in the LOWL 
set - have also been introduced. 



4.3 The "best" mode set 

Having ascertained that the BiSON and LOWL data sets 
were compatible, we selected a "best" low-degree mode set 
for the final inversion results from the four contemporaneous 
sets (sets 1-3, 5 of Table 1). The mode selection criteria 
were as follows. If a mode was present in one set only it was 
selected automatically. This was true for most of the modes 
in the 32-month set. When modes were present in both the 
BiSON and LOWL sets, the determination with the lower 
error was selected. When the errors were essentially equal, 
the average frequency was used, together with the error from 
just one of the sets, i.e., we took a conservative approach and 
did not reduce the error when averaging, on the grounds that 
the errors in the two sets were not independent (there being 
a correlated contribution due to solar noise). When faced 
with a choice between the BiSON-8 and BiSON-G data, the 
BiSON-G data were used if the fit to the BiSON-8 mode 
was obviously poor. For I = 4 and above, the LOWL data 
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Figure 5. The low-degree part of the l — i/ diagram of the "best" 
set. The filled triangles are LOWL modes; circles are BiSON-8 
modes; asterisks are modes from the 32-month set; filled squares 
the BiSON-G modes; and empty triangles averages of LOWL and 
BiSON-8 modes. The error bars represent lOOOcr errors. 
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Figure 6. The sound-speed inversion results of the "best" set, 
in the sense (Sun) - (Model). The lower panel (b) concentrates 
on the results in the core. For comparison, the la bounds on the 
result of inverting just the LOWL set are shown by the dotted 
line. The plot is on the same scale as Fig. 4. 

were used. The I < 3 part of the l-v diagram for the Best 
Set is shown in Fig. 5. 

The sound-speed inversion results of the Best Set are 
shown in Fig. 6. The main change in the inversion, over 
the LOWL-only calculation, is that the propagated errors 
in the core inversion have been reduced substantially, i.e., 
by more than a factor of two at many radii. The propagated 
errors in the inversion for the LOWL-only set and Best Set 
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Figure 7. A comparison of the propagated errors and resolution 
of the inversion results of the LOWL and "best" mode sets. The 
upper panel shows the la errors on the solution at different target 
radii. The dashed line is the error on the inversion results for the 
LOWL set, while the continuous line is that for the "best" set. In 
the lower panel, the thick continuous line is the width (in units of 
the solar radius Rq) of the averaging kernels for the best set. The 
thick dotted, thick dashed and thick dot-dashed line are the differ- 
ences between the first, second and third quartile points and the 
target radius for the best-set inversion results, all in units of Rq . 
The thin line shows the corresponding quantities for the inversion 
with the LOWL data. The width of the averaging kernel, defined 
as the distance between the first and third quartile points. Note 
that the averaging kernels become symmetric at around O.IRq. 

are compared in Fig. 7. Note that the influence of the low- 
degree modes on the errors is restricted to radii less than 
0.35 R©. A few averaging kernels have been shown in Fig. 8. 

As can be seen from Fig. 7, the data combination also 
results in a slight improvement in resolution. This can be 
clearly seen by inverting the frequency differences between 
the reference and the test model - using only the available 
modes in each set and weighting by the observed errors. The 
results of the inversions with the Best Set and LOWL-only 
set are shown in Fig. 9. The sound-speed difference between 
the two models increases towards the core. The inverted dif- 
ference is an average of the exact difference over a flnite ra- 
dius - consequently, unless the averaging is performed over 
a narrow region, the inverted difference cannot match the 
exact difference. Thus a higher-resolution inversion will give 
a better match to the exact differences. While the results 
are, as expected, identical in the outer layers, the inversion 
results in the core are not. The results of the Best Set are 
marginally closer to the exact model differences than are 
those of the LOWL set alone. 

We have also carried out inversions to determine the 
relative density difference between the Sun and the refer- 
ence model, using the pair of variables {p,Y). The results, 
using both the LOWL-only set and Best Set, are shown in 
Fig. 10. Unlike the sound-speed inversion results, we find 
that the density results for the two data sets differ at all 
radii, even though the sets differ only in low-degree modes 
which predominantly carry information about the core. Also, 
the differences in 5p/p are substantially larger than the es- 
timated random errors in the results. Probably the explana- 
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Figure 8. The averaging kernels at target radii of (a) 0.063 Rq 
and (b) 0.51 Rq. The continuous lines indicate the best-set in- 
version and the dotted lines the LOWL inversion. Marked on the 
averaging kernels are the three quartile points. The asymmetry 
between the positions of the quartile points indicated that the av- 
eraging kernels have some structure away from the target radius. 

tion for these properties is that our (p, y) inversion is sensi- 
tive to discrepancies between the reference model's equation 
of state, which was used in deriving the (p, Y) kernels, and 
the true equation of state of the solar plasma (cf. Basu & 
Christensen-Dalsgaard 1997). Note that in this context these 
discrepancies can also include the effect on Fi of differences 
in the abundance of heavy elements. In additional tests with 
artificial data we have found comparable effects of using the 
two different mode sets if the test model was based on the 
so-called MHD equation of state (e.g. Mihalas, Dappen & 
Hummer 1988) while the reference model used the OPAL 
equation of state. It seems credible that the intrinsic error 
in the OPAL equation of state, relative to the Sun, is of 
a magnitude similar to the difference between OPAL and 
MHD. We also note that the mass-conservation constraint. 



&P 2, 
p — r dr 

P 







(5) 



implies a strong correlation between the results of density 
inversions at all radii. Thus even if the change in mode set 
has its primary effect in the core, the inferred density dif- 
ference is likely to change at all radii, as observed. Such 
sensitivity to the equation of state can be greatly reduced 
by using (p, Fi) or by explicitly seeking to suppress in the 
(p, Y) inversion any contributions from intrinsic differences 
in Fi (Basu & Christensen-Dalsgaard 1997), although this 
results in larger random errors in the solution for reasons we 
have already discussed in Section 1. 

As in the case of the sound-speed inversion, the errors 
in the density inversion are reduced and the resolution im- 
proved by using the Best Set, compared with inversion of 
the LOWL-only set. Also, the tests on artificial data sug- 
gest that the effects of errors in the equation of state are 



© 1994 RAS, MNRAS 000, 000-000 



8 S. Basu et al. 



0.004 



0.002 




^«n*-^^- 



-o.ooa 



0.2 



0.4 



0.6 



0.8 



0.015 - 



0.01 



a 0.005 - 



-0.005 



0.01 















0.2 



0.4 



0.6 



0.8 



r/R„ 



r/R„ 



Figure 9. The continuous line stiows tlie exact sound-speed dif- 
ference between tlie test (4.52 Gyr) and reference (4.6 Gyr) mod- 
els, in the sense (younger) minus (older). The symbols show the 
results of the inversion to determine the difference between mod- 
els. The circles are inversion results using the Best Set modes 
and the triangles those using the LOWL modes. Note that near 
the core, the results using the Best Set are closer to the exact 
difference. This indicates better resolution. 

less severe for the Best Set than for LOWL set, probably as 
a result of the stronger constraints on the core provided by 
the former set. 



5 DISCUSSION AND CONCLUSIONS 

We have shown that solar-oscillation data from different 
sources can be successfully combined if the observations 
are contemporaneous. Although a homogeneous dataset is 
preferable, different instruments are best suited to observing 
different ranges of modes. In particular, observations made 
in integrated sunlight - which are sensitive only to modes of 
the lowest degrees - can be combined with resolved obser- 
vations to give a reasonably large selection of modes. 

Using contemporaneous data means that corrections for 
solar cycle effects - which are still very uncertain - are not 
required. The results obtained from such a combined set are 
consequently more reliable. The combination of data from 
the LOWL instrument and the BiSON network shows that 
the inversion errors in the core can be reduced by a up to a 
factor of two in the part of the Sun whose structure is the 
most difficult to determine. The resolution of the inversions 
is also improved to a certain extent. 

We find from Fig. 6 that the sound speed c of a solar 
model constructed with the most up-do-date physics avail- 
able is very close to that of the Sun. The difference in (? 
is at most 0.5 per cent and is considerably less in most 
regions. However, this small discrepancy is quite substan- 
tial compared with the errors in the inversion, as inferred 
by propagating the data errors. The sound-speed difference 



Figure 10. The relative density difference between the Sun and 
the reference model, in the sense (Sun) - (Model), as inferred 
by the inversion of the LOWL (triangles) and the "best" sets 
(circles). 

between the Sun and the model has three very noticeable 
features — a large positive bump at roughly O.TR©, a dip 
around O.2R0 and a rise in the core. The first feature is just 
beneath the base of the solar convection zone, which lies at 
roughly O.7R0 (Christensen-Dalsgaard, Gough & Thomp- 
son 1991; Basu & Antia 1997). The sound speed is given 
by (? — Tip/p oc TiT/fi in the excellent approximation of 
an ideal gas. Thus the presence of the first feature could be 
explained by helium settling being too strong in the model 
in this region. The increase in the helium abundance in- 
creases the mean molecular weight of the material at the 
convection-zone base, and hence decreases the sound speed 
in the model. If the gravitational settling has been overes- 
timated, the reason could be that in the Sun the settling is 
inhibited by some mixing in this region (Gough et al. 1996), 
a speculation supported by secondary inversions for the he- 
lium profile in the Sun (Antia & Chitre 1997). A mismatch 
in the depth of the convection zones in the model and Sun 
(Basu & Antia 1997) could also contribute to the feature. 

Probing the structure of the energy-generating core pro- 
vides a particularly important observational test of the the- 
ory of stellar evolution. The negative region and the sub- 
sequent rise in the sound-speed difference implies that the 
variation of the sound speed in this region of the Sun is 
flatter than in the reference model. In this region the abun- 
dance varies with position not only because of settling but 
also due to the burning of hydrogen into helium over the 
Sun's lifetime. The inversion results could be explained if 
some mixing has transported helium outward from the in- 
ner core, leading to the Sun's sound speed being increased 
in the inner core and decreased around 0.2R, relative to the 
model (cf. Gough et al. 1996). 

The density differences between the Sun and the refer- 
ence model, illustrated in Fig. 10, are also of modest size, 
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although substantial compared with the estimated errors. 
Overall, the density in the Sun is higher than that in the 
model, except in the core - the Sun is therefore less centrally 
condensed than the model. The density difference shows a 
jump at the base of the convection zone, corresponding to 
the sharp feature in the sound-speed difference at this lo- 
cation. Density differences within the convection zone are 
more sensitive to model parameters than the sound speed. 
Such discrepancies could be caused by differences in con- 
vection zone depths and elemental abundances. The results 
for density show a considerable, and somewhat worrying, 
sensitivity to the choice of mode set. We have argued that 
might plausibly be caused by errors in the equation of state 
assumed in the calculation of the reference model. Such ef- 
fects can in fact be suppressed in the inversion, although at 
the expense of a substantial increase in the random errors 
(cf. Basu & Christensen-Dalsgaard 1997). On a more pos- 
itive note, the sensitivity hints at the potential to extract 
even quite subtle information about the equation of state 
from helioseismic inversions. 

The positive sound-speed difference in the core and the 
smaller central condensation might indicate that the Sun is 
younger than the assumed age of 4.6 Gyr. Indeed, estimates 
based on meteoritic ages (e.g. Guenther 1989; Wasserburg, 
in Bahcall & Pinsonneault 1995) indicate that the age is 
close to 4.52 Gyr, as assumed in the test model. Comparison 
of the exact difference in Fig. 9 with the inversion in Fig. 6 
shows that using a younger reference model would tend to 
reduce the difference between the Sun and the model very 
near the centre, but it would increase it elsewhere. Similarly, 
Basu et al. (1996a) found that although the density differ- 
ence between the Sun and a younger model is lower in the 
core, the differences elsewhere increase. Alternatively, some 
degree of mixing in the core could resolve both the problem 
of the larger sound speed and the lower density in the solar 
core. 
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